CRACK INITIATION IN ELASTIC BODIES 



ANTONIN CHAMBOLLE, ALESSANDRO GIACOMINI, AND MARCELLO PONSIGLIONE 

Abstract. In this paper we study the crack initiation in a hyper-elastic body governed by 
a Griffith's type energy. We prove that, during a load process through a time dependent 
boundary datum of the type t — > tg(x) and in absence of strong singularities (this is 
the case of homogeneous isotropic materials) the crack initiation is brutal, i.e., a big 
crack appears after a positive time ti > 0. On the contrary, in presence of a point x 
of strong singularity, a crack will depart from x at the initial time of loading and with 
zero velocity. We prove these facts (largely expected by the experts of material science) 
for admissible cracks belonging to the large class of closed one dimensional sets with a 
finite number of connected components. 

The main tool we employ to address the problem is a local minimality result for the 
functional 

£{u,Y) := I f(x,Vv)dx + kH 1 (T), 
Jn 

where Q C R 2 , k > and / is a suitable Caratheodory function. We prove that if the 
uncracked configuration u of Q. relative to a boundary displacement ip has uniformly 
weak singularities, then configurations (ur,r) with 7i 1 (F) small enough are such that 
£{u,0) < £{u r ,T). 

Keywords : free discontinuity problems, energy minimization, crack initiation, varia- 
tional models. 
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1. Introduction 

Griffith's criterion for crack propagation in hyper-elastic bodies asserts that, during a load 
process, a crack T can grow only if the energy dissipated to enlarge the crack, which is basically 
assumed to be proportional to the area of the cracked surface, is balanced by the corresponding 
release of bulk energy. According to Griffith's theory, if Q represents a two dimensional hyper- 
elastic body, tjj is a boundary datum and T is a curve in Q parametrized by arc length, then the 
crack r(/o) is in equilibrium if 

nr] ,, M r W(u(Iq)) - W(u(l + Q) 
(1.1) k(lo) ■= hmsup : < k, 
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where u(Iq) and u(Iq + 1) are the displacements associated to ip and to the cracks T(l ) and T(l n + l) 
respectively, W is the bulk energy functional and k is the toughness of the material. A quasistatic 
crack evolution is determined by an increasing function t — » l(t) satisfying the Griffith's criterion 
for crack propagation, which asserts that for every t we have k(l(t)) < k and 

(k - k(i(t)))i(t) = o, 

i.e., T(l(t)) propagates if and only if holds with equality. 

The aim of this paper is to discuss this criterion in the case of crack initiation, i.e., when there 
is not a pre-existing crack in the body (7o = 0). A fundamental role in the problem is played 
by the singularities of the body, namely the behavior of the elastic energy concentration of the 
deformation. Experiments show that small cracks usually appear near sufficiently strong singular 
points of the the body, whose position arc essentially determined by its inhomogcneitics. If the 
singularities of the body are sufficiently weak (for instance this is the case of homogeneous isotropic 
materials), a lot of results in the literature of material science show that the derivative in for 
lo = is equal to zero. The conclusion is that Griffith's criterion is not adequate to predict crack 
initiation (and, as a consequence, a crack evolution originating from an uncracked configuration). 
These results require that the path of the crack is sufficiently regular (a line or a smooth curve). 
In this paper we prove that the same conclusion holds in the class of all one dimensional closed 
sets with a finite number of connected components. More precisely we prove that the limit in (|l.lfl 
is zero if is any family of closed sets with length less than I and with at most m connected 
components, with m independent of I. In particular we do not prescribe the path nor the shape 
of the cracks. 

Although it is more general, our study is in part motivated by the variational model for qua- 
sistatic crack propagation proposed by Francfort and Marigo in [T^J. The main features of this 
model are that the path of the increasing crack T(t) is not preassigned, the class of admissible 
cracks is given by all sets with finite length, and the growth is not assumed to be progressive, 
namely the length of the crack is not assumed to be continuous in time. The classical Griffith's 
equilibrium condition for the configuration (u(t), T(t)) is replaced by a unilateral minimality prop- 
erty and an energy balance condition. The unilateral minimality property states that, during the 
crack evolution, the total energy is minimal among all configurations with larger cracks, namely 

(1.2) W(u(t)) + fcW x (r(t)) < W(v) + kn\H), 

for every crack H larger than T(t) and for every deformation v admissible for the boundary datum 
ip and for H . (Here Ti 1 — the 1-dimensional Hausdorff measure — is a suitable generalization of 
the length). The energy balance condition states that the energy of the system evolves in relation 
with the power of external loads in such a way that no dissipation occurs. The authors claim 
that their model improves the understanding of the crack initiation with respect to the classical 
Griffith's criterion: as a matter of fact, in contrast with Griffith's model, it admits brutal crack 
initiation, i.e., evolutions T(t) of the type 

r(i) = for every t < U 

and 

inf ttVlY*)) > 0, 

t>ti 

where ti is referred to as time initiation of the crack. In this paper we prove that, within the class 
of cracks which are closed and with at most m connected components, crack initiation is always 
brutal whenever the elastic displacement presents sufficiently weak singularities. More precisely 
we show that (|1.2fl can be violated only by cracks whose length is greater than a critical quantity 
I*, depending on the boundary datum and on the physical properties of the material. On the 
contrary, in presence of a point x of strong singularity the crack initiation is progressive: a crack 
departs from x at the initial time of loading and with zero velocity. These facts were proved by 
Francfort and Marigo in Proposition 4.19], under the assumption that the path of the crack 
is given a priori by a finite number of fixed curves which can be parametrized by arc length. This 
is not the case in our larger class of admissible cracks. 
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The main tool we employ to address the problem of crack initiation is a local minimality result 
for the functional 

(1.3) / f(x,Vv)dx + kn 1 (T), 

Jn 

where fl is a bounded Lipschitz open set in R 2 , k > 0, and / : fl x R 2 — > R is a Caratheodory 
function strictly convex and C 1 in the second variable, satisfying standard p-growth estimates 
with p > 1, and such that f(x,0) = 0. The functional (|1.3|) is a variant of a functional which 
has first appeared in the theory of image segmentation, in a celebrated paper by Mumford and 
Shah 20 . The set T belongs to the class 

(1.4) /C m (fl) :— {r C fl : T has at most m connected components and Ti. l (T) < +00} 

and the function v belongs to the Sobolev space W 1,P (U\T) and satisfies the boundary condition 

(1.5) v = 4> ond D n\T, 

where <9ofl C <9fl is open in the relative topology, and ip is (the trace of a function) in W 1 ' P (SY) fl 
L°°(fl). 

Let Mr be a minimum energy displacement relative to ip and T, i.e., let ur be a minimizer for 

(1.6) minjy f(x,Vv)dx : u £ W^ p (fl\r),v = ip on d D fl 

We denote by u the elastic configuration of fl relative to the boundary datum ip, i.e., a solution 
of l|1.6fl with r = 0, and we assume that u admits uniformly weak singularities in fl, i.e., 

(1-7) HV«||£, (Br nn) < Cr a 

for some constants a > 1 and C > and for every ball B r with radius r. Condition i|1.7fl means 
that the bulk energy of the elastic configuration u in a ball B r {x) is negligible with respect to the 
length of dB r (x) as r goes to zero, uniformly in 1 e fl. 

Our main result is the following Theorem, which establishes that under the previous assumptions 
small cracks are not energetically convenient for the functional (|1.3|) . 

Theorem 1.1. Assume that u admits only uniformly weak singularities in Q. Then there exists a 
critical length I* depending only on fl, /, k, ip and m such that for all T G /C m (A) with H 1 (r) < I* 
we have 

(1.8) / f(x,Vu)dx< I /(s,Vu r )rfx + fcW 1 (r). 



We observe (see Remark 13. 311 that this statement is equivalent to the local minimality of u 
in l|1.3fl . in the L 1 topology. 

Let us briefly comment the assumption about the singularities of u in Theorem 11.11 The 
minimality result is false if the elastic solution u has strong singularities, namely if there exists 
x £ fl such that 

(1.9) limsup- / \Wu\ p dx = +00. 

r^o r JB r (x)nn 

In fact condition ljl.9|l ensures that it is energetically convenient to create a small crack T := 
dB r (x) around x: the surface energy needed to create such a crack is proportional to r, while the 
corresponding release of bulk energy is by l|1.9fl bigger than r if r is small enough. 

The critical case when the right hand-side of (|1.9(l is a constant < C < 00 corresponds to 
the singularity appearing around the tip of the crack (see |14j). In this case the celebrated Irwin's 
formula states that the release of bulk energy per unit length along rectilinear increments of the 
crack is equal to the so called mode III stress intensity factor Km, which is proportional to C . In 
our class of cracks /C m (fl) we have that if C is small enough, then the release of bulk energy per 
unit length is less than k, and therefore our minimality result still holds, while it is false if C is 
too large. We can not fill the gap, and therefore we do not achieve a sharp Irwin type formula in 
our class of cracks. 
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In order to prove Theorem 11.11 we have to compare the asymptotic behavior of the release of 
bulk energy 

(1.10) f[f(x,Vu)-f(x,Vu r )]dx 

Jn 

with H x (r) when « x (r) -> 0. In the literature there are many results in this direction considering 
particular sequences of infinitesimal cracks r n , for instance when T n is the rescaled version of a 
fixed smooth curve T. An intuitive strategy to estimate (jl . is to compute how much energy is 
required in order to make Mr a good competitor for the minimum problem p. 6(1 without cracks, 
namely how much energy is required to heal the crack T. 

This seems difficult for a generic crack in K. m (Q). So our strategy is to operate on the stress 
(7 := df(x,S7u) of the elastic solution. More precisely we prove the following key estimate (see 

(1.11) f [f(x,Vu)-f(x,Vu r )]dx< f [T-a}-[dr(x,T)-df*(x,<r)]dx, 
Jn Jn 

for all vector fields t S I«(fl;l JV ) (q = p' := -^-) such that 

(1.12) / r-V« = 0forall«e A{T). 

Jn 

Here /* is the convex conjugate of /, and A(T) is defined as 

A(T) := {v e W hp (n \ T) : v = on d D n \ T}. 

If r is sufficiently regular, condition l|1.12fl means that the vector field r has zero divergence outside 
r, and t{x) is tangent to T for every x £ T. 

The proof of Theorem 11.11 relies on the construction of a vector field r satisfying l|1.12fl and 
such that 

(1.13) f [t - a] ■ [df*(x,T) - df*(x,cr)}dx < kH 1 ^). 
Jn 

We construct r modifying the stress a which is divergence-free (as a consequence of Euler equation 
of problem l|1.6|l ). but not tangent to T. First, we consider a neighborhood U of T, and a cut-off 
function (p such that (p = on U. Then ipa is null near T, and in particular it is tangent to T. 
Then we construct a vector field r\ in such a way that 77 = on U and 

divr] — — div(^cr). 

We get that r = ipa + r\ is an admissible vector fields for inequality (|1.11|) . Using the fact that 
r G /C m (f2) and that u has uniformly weak singularities, it is possible to choose U, (p and rj in 
such a way that inequality i|1.13[l holds. 

It turns out that the constraint l|1.12fl can be handled in an easier way than the constraint of 
being a gradient, and this is the reason why we work with the stress a instead of the strain Vu. 

Estimate (|l.lll) actually holds true in any dimension, and it turns out that our arguments work 
in any dimension provided that the crack T belongs to the class 

(1.14) /C C (H) := {r C II : T is closed and diam(r) < CH N ^ 1 (T)}, 

where C is a fixed constant. This is certainly true in dimension two for the cracks in /Ci(f2) (and 
in /C m (0) up to a localization argument). However, it also shows that the local minimality result 
of Theorem 11.11 remains valid in higher dimension, within the class IC C (fl) of cracks that are not 
needle-like: see Remark 13. 121 

The minimality result holds also in the case of planar linearized elasticity, with a density of bulk 
energy involving the symmetrized gradient. This is considered in Section 0] — while the "simpler" 
case of 2D vectorial nonlinear elasticity is addressed in Remark 13.111 

Finally in Appendix lAl we extend the local minimality result of Theorem II. II to the larger class 
of all 1-dimcnsional rectifiable sets. This seems to be the most general class: however, there is 
a price to pay in order to handle such admissible cracks. First, we need to assume that Vu is 
bounded and regular up to the boundary. In particular we are not able to treat the case in which 
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u has weak singularities. Then, the method we employ is based on the maximum principle, which 
allows to estimate the local opening of a crack with the global energy in a small ball surrounding 
the crack. It is therefore strictly scalar, and bidimensional. 

The paper is organized as follows. In Section [21 we establish the main inequality In 
Section |3| we prove the local minimality result in dimension 2 and its extensions to vector- valued 
displacements and to the iV-dimensional case within the class of cracks given by l|1.14fl . Sec- 
tion 0] addresses the case of planar two-dimensional elasticity. The problem of crack initiation 
in quasistatic evolutions is addressed in Sectional while the two-dimensional SBV-case without 
singularities is treated in Appendix In Appendix [51 we show how to obtain some uniform 
Poincare and Poincare-Korn type inequalities, used during the proofs of our main results in the 
construction of the competitor stress field r. 

2. The dual problem and the main estimate 

Let fi be a bounded connected Lipschitz open set in R , let 9dQ C dft be open in the relative 
topology, and let On^I := dfl \ <9of2. Let / : x R N — > R be a Caratheodory function such that 

(2.1) £ — > /(&,£) is strictly convex and C 1 for a.e. xe!!, 



(2.2) f(x,0) = for a.e. x G O, 
and such that for a.e. x G Q and for all £ G 

(2.3) a\^<f{x^)<m\ P + ^), 

where a, (3 > and 1 < p < +oo. 

Given ip G W 1,p (fi) fl L°°(f2) and T closed set contained in fl, let us consider the minimization 
problem 

(2.4) min j^J f(x, Vu) dx : u G W hp (Cl \ T), u = ip on d D n \ r| . 

In view of and since ip G W 1,p (Cl) fl L°°(17), problem l|2.4|l is well posed. 

Let us denote by ur G W /1 ' p (fi \T) n £°°(£1 \T) a minimizer of (j2.4jl . Clearly Vitr is uniquely 
determined, while ur is determined up to a constant on each connected component of 51 \ L which 
does not touch du^i- 

We denote by u G W 1,p (tt) n L°°(fi) the solution of (|2.4H corresponding to L = 0, and we 
refer to u as the elastic solution. In the case T is sufficiently regular, the Euler Lagrange equation 
satisfied by u-p is 

!div9 c /(x,Vur) = on fl\T, 
u = ip on d D fl \ T, 

Vw r ) ■ n = on d N Q U <9I\ 

where n denotes the normal vector to d/vfl LIT, In the sequel, we will write df(x,t;) for d^f(x,^). 
Let us set 

(2.6) A(T) := {v G W 1 *^ \ T) : « = on \ T}. 

Let us denote by /* the convex conjugate of / with respect to the second variable defined by 

f(x,{) : sup{C-£ /U.O : £GR W }. 

We refer the reader to [21] for the main properties of the conjugate function /*. Notice that /* is 
of class C 1 since / is strictly convex. The main result of this section is the following. 

Theorem 2.1. Let T be a closed subset ofQ, and let a := df(x,\7u) be the stress associated to 
the elastic configuration u. Then we have 

(2.7) / [f(x, Vu) - f(x, Vu T )} dx < [ {t ~ a] ■ [df*(x, r) - df*(x, a)] dx 
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for all t G Li{fl; R N ) (q = p' := ^) such that 

(2.8) / t -Vvdx = for all v G A(T). 

Jn 

Proof. For all 77 £ LP{VL; R n ) let us set 

(2.9) $(??):= min / f(x,Vw + rf)dx, 

weu+A{T) J n 

where -4(r) is defined in (|2.6(l . Then the convex conjugate of $ defined on L q (fl;M. N ) has the 
form 

(2.10) $*(r) := sup / [r • 77 - /(at, Vu; + 77)] dx 

™,»? in 

= sup / [r • (77 + Vw) — /(x, Vic + 77) — t ■ Vw] dx 

— / [f*{ x i t) — t ■ Vu] da; + sup / r • Vv dx. 
in «eX(r) Jn 

We conclude that 

f f\f*(x,T) -T- Vuldx if t satisfiesESJ) 

(2.11) **(t) = ^ in V ; J ' 

l+oo otherwise. 
Notice that ^(O) = $**(0) because $ is weakly lower semicontinuous and 

$(0) = / Vu r )dx < +00. 
is) 

(In fact, $ is locally finite, hence locally Lipschitz). Therefore we obtain 

- / /foVu r )dz = -$(0) = -$**(0)=rmn$*(r), 

so that by (|2.11|) we deduce 

— / /(x, Vu r ) dx — min ^ / [f*(x, t) — t ■ Vu] dx : r satisfies l|2.8l) > . 
Jn T Un J 

For all t satisfying JjOJ, we get 

[/(x, Vu) - f[x, Vu T )} dx < I [/(x, Vu) + f*(x, t) - t ■ Vu] dx. 
2 Jn 

Let cr(x) := df(x, Vu(x)) be the stress of the elastic solution u. Since for a.e. x € fl 

f*(x,T(x)) < r(x,a(x))+dr(x,r(x)) ■ (r(x) - a(x)), 

and 

f{x, Vu[x)) + f*(x, a{x)) = Vu{x) • a(x), 

and since Vu(x) — df*(x,a(x)), we finally obtain our main estimate l|2.7|) . so that the proof is 
concluded. □ 

Remark 2.2. If T is sufficiently regular, condition (|2.8|) means that r has zero divergence outside 
r, and t(x) is tangent to T for every x G dT. 
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3. The minimality result in anti-plane elasticity 

In this section we prove that under some assumptions on the elastic configuration u, small 
cracks are not convenient for the total energy 

(3.1) / f(x,Vv)dx + kH 1 (T), 

Jn 

where / is a Caratheodory function satisfying conditions (|2.1|l . (|2.2|) and (|2.3(1 . and k > 0. 

Let us consider ft bounded connected Lipschitz open subset of R 2 , and let dofi C d£l be 
open in the relative topology and such that <9/v^ := dfl \ 9of2 has a finite number of connected 
components. 

Let to be a positive integer. The class of admissible cracks is given by 

(3.2) /C m (f2) := {L C Q : T has at most to connected components and 7i 1 (L) < +oo}. 
Given ip € W 1 ^^) n L°°(f2) and T € /C TO (fi), the displacement 

u T e w 1,p (n\r) 

associated to L and ip is given by problem l|2.4|l . We denote with u the solution of l|2.4ll relative 
to r = 0, and we refer to u as the elastic solution. 

The basic assumption on the elastic configuration u involves the behavior of the energy con- 
centration of the stress. We require that u has uniformly weak singularities in f2, in the sense of 
the following definition. 

Definition 3.1. We say that u £ W ,p (f2) has only uniformly weak singularities in A, open subset 
of il, if there exist constants 1 < a < 2 and C > such that for every x S A and for every r small 



(3.3) / |Vu| p dx < Cr a . 

As mentioned in the Introduction, the condition of uniformly weak singularities means that the 
bulk energy of the elastic configuration u in a ball B r (x) is an infinitesimal of higher order than 
the length of dB r (x) as r goes to zero, uniformly with respect to x G A. 

The main result of this section is the following theorem. 

Theorem 3.2. Let the elastic solution u of problem ()2.4() have only uniformly weak singularities 
in n according to Definition ^. 11 Then there exists a critical length I* > depending only on fi ; 
m, f , k and ip such that for all T £ IC m (il) with 7i 1 (r) < I* we have 



(3.4) 



[ f{x,X7u)dx< [ f[x,Vu T )dx + kH 1 {T), 
Jn Jn 



where ur is a minimum of Q2.4JI ■ 

Remark 3.3. (Local minimality in the L 1 -topology) The minimality condition (|3.4() implies 
that the elastic solution u is a local minimum for the total energy l|3.1|l with respect to the L 1 - 
topology. More precisely, for every sequence (T^heN m JC m (fl) and for every Uh & W 1,p (fl \ T/J 
with Uh — > u strongly in L 1 (il), for h large enough we have 

(3.5) / f{x,x/u)dx< I f{x.yu h )dx + kn\r h ). 

Jn Jn 

In fact, it is not restrictive to assume that (Vu^h^N is bounded in L p (fi, R 2 ) and that TC 1 (T^) < C . 
By Ambrosio's lower semicontinuity theorem [2] we have 

/ f(x, Vu) dx < liminf / f(x,Vui l )dx. 
Jn h^r+co J n 

If the sequence (7i 1 (r/ l ))^ e N is not infinitesimal, then (|3.5(l clearly holds. If 7i 1 (r^) — > 0, we have 
W 1 (r/i) < I* for h large enough, and hence l|3.5l) follows from l|3.4ll . 

Actually, using Ambrosio's compactness theorem [2], one can show that the local minimality in 
the I^-topology and the minimality result of Theorem 13.21 are equivalent. 



8 



A. CHAMBOLLE, A. GIACOMINI, AND M. PONSIGLIONE 



In order to prove Theorem 13 . 21 we will use the main estimate given by Theorem 12. II Our aim 
is to construct a vector field r S L q (il;M. 2 ) which is an admissible competitor in l|2.7|l and which 
shows that the difference between the bulk energies of u and ur is smaller than kTi}{T). In order 
to do so, we need some preliminary lemmas. 

Lemma 3.4. Let x € Q, and let r > be such that B2 r {x) C Q. Let ip be a smooth function 
with 0<ip<l, ip = 0on B r (x), cp = 1 outside B2 r {x) and HVtpHoo < -. Then there exists 
f] £ L q ((B 2r {x) \B r (x));R 2 ) with q := p' = p/(jp- 1) such that 

J div 77 = —div((pa) on B2 r (x) \ B r (x), 
I 77 • n = on d(B 2r {x) \ B r (x)) 



and 

' ' \a\ q dx 



(3.7) 



(3.6) f \r,\ q dx < C [ 

J B 2r {x)\B r {x) J B 2r (x)\B r (x) 

where a := df(x,Vu) is the stress of the elastic solution u, n is the outer normal to d(B2 r (x) \ 
B r {x)), and C is a constant independent of r. 

Proof. Let us set 77 := |Vu| p ~ 2 Vt>, where v g W 1 ' p (B2 r {x) \ B r (x)) satisfies the equation 

jdiv(\\7v\P- 2 \7v) = -div((pa) in B 2r (x) \ B r (x), 
\ \X7v\p- 2 X7v ■ n = on d(B 2r {x) \ B r (x)). 

Notice that the equation is well posed because 

(3.8) -div(^o-) = -Vip ■ a e L q {n) 

and 

/ div((pcr) dx = / a-ndTL 1 = Q 

J B 2r (x)\B r (x) JdB 2r {x) 

since a is divergence free in Q. 

It remains to prove inequality l|3.6|l . To this aim, note that we can always assume that v has 
zero mean value, so that by Poincare inequality and by a rescaling argument, we have that there 
exists C > independent on r such that 



(3.9) / \v\ v dx < Cr p / \Vv\ p dx. 

J B 2r (x)\B r (x) J B 2r {x)\B r (x) 

Recalling that HV^Hoo < 2/r and taking into account l|3.8|) and l|3.9|l we get 

/ \\7v\ p dx= / ^(fi-vvdx <\\Vtp-a\\ Lq{ B 2r (x)\B r (x))\\v\\Lp(Bo r (x)\B r (x)) 

J B 2r (x)\B r (x) J B 2r (x)\B r {x) 

< -\\< T \\L«(B 2r {x)\B r (x)M 2 )C 1 ^ Pr \\V" v \\Li>(B 2r (x)\B r (x),R 2 ) 

= C 1/p ||fT|| Lg(B2r ( I )\ Br ( a; ), R 2)||\77;|| L p( B2r ( :l; )\ Br ( a; ) iR 2) 

so that 

f \rj\ q dx= [ \Vv\ p dx < C q/p [ \a\ q dx, 

J B 2r (x)\B r (x) JB 2r (x)\B r (x) J B 2r (x)\B r (x) 

and this concludes the proof. □ 

Now we need to construct a suitable field 77 also around points x on the boundary of f2. Since 
n is Lipschitz, for every x S dfl we can find an orthogonal coordinate system (a^,^) with origin 
at x, £i,£2 > and a Lipschitz function g : [— £i,£i] — > [—£2,^2] such that setting 

(3.10) R r (x) :={(xi,4) : \x[\ < re u \x' 2 \ < re 2 } 
we have for r small enough 

(3.11) ft n R r (x) = {(x'i, x' 2 ) G R r (x) : 4>fif(a;i)}. 
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Notice moreover that the Lipschitz constant C g of the function g is determined only by f2, and 
that we can assume C s e\ < e^. Let us set 

(3.12) A r (x) := (R 2r (x)\R r (x))nn. 

Lemma 3.5. Let r > be small enough (so that (|3.11|) holds everywhere on dfl). Let x G dfl 
such that one of the following three situations holds: 

(1) x G d D £l and A r (x) H dfl C d D Cl; 

(2) x G int(9jvf2) and A r (x) D dfl C diyQ (where int(-) indicates the interior relative to dQ); 

(3) x G dfl \ {d D fl U int{d N n)). 

Let p be a smooth function with < p < 1, p — on R r (x), tp = 1 outside R2r{x) and WV(p\\oc < r- 
Then there exists n G L q (A r (x); M 2 ) with q := p' = pj (p — 1) such that 

{div 77 = — div(ipo-) in A r (x) 
rj-n = ondA r (x)D(nud N n) 

and 

[ \v\ q dx < C I \a\ q dx, 

JA r {x) JA r (x) 

where a := df(x, Vu) is the stress of the elastic solution u, n is the outer normal to dA r (x) n (f2U 
djyQ), and C is a constant depending only on fL 



We observe that only the points in d£l at distance less than rysf + ef to dCl\ ((9£>f2Uint(<9jv^)) 
might not fall into one of first two cases. 

Proof. In all the three cases of the lemma, we will use the fact that the Poincare inequality holds in 
W 1,p (A r (x)) with a constant that rescales as r, i.e., there exists a positive constant C independent 
of r such that 



(3.13) / \v\ p dx < Cr p / \Vv\ p dx 

JA r {x) J A r (x) 

for all v G W 1 ^{A r {x)) with 

v = or v = on OdO, Pi dA r {x). 

A r (x) 

This can be seen rescaling A r (x) with the transformation T r (x',y') = (fr, ^r), and using Propo- 
sition ^Q] (see Appendix lB|l in the domains {T r (A r (x))}: it shows that the Poincare inequality 
holds in T r (A r (x)) with a constant that is independent of r, and by rescaling we deduce that 
(|3~T3|) holds. 

In case (1), we can consider 77 := |Vw| p_2 Vw, with v G W 1 ' p (A r (x)) satisfying the equation 



(3.14) 



' div(\Vv\ p - 2 Vv) = -div(<pa) in A r (x), 
v = on dA r {x) n d D il, 

]Vv\ p - 2 Vv-n = Q on dA r (x)nn. 



Notice that — div^er) = — Vp ■ a (a is divergence free) and || V<£>||oo < 7- Taking into account the 
Poincare inequality i|3.13[) we get 



\Vv\ p dx= / Vp ■ avdx < \\Vp ■ o-\\ LqiAr{x ))\\v\\ LP(AAx)) 

A r (x) JA r (x) 



< -\\ a \\L-i(A r (x).R 2 )C 1/Pr \\^ v \\LP(A T (x)M 2 ) = C 1/P \\< 7 \\L''{A r (x)M 2 )\\^ v \\LP(A r (x),R 2 ) 

[ \r 1 \ q dx = [ \Vv\ p dx < C q/p [ \a\ q dx. 

JAJx) J AJx) JAt(x) 
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Cases (2) and (3) can be treated as case (1) considering ?y := \Vv\ p 2 Vi> with v € W 1,p (A r (x)) 
denned by the equations 

Uiv(\Vv\P-^v) =-dw( W ) \xvAr(x), 
1 ' ' \\Vv\ p ~ 2 Vv ■ n = ondA r {x), 

and 

'div(|Vv| p - 2 Vu) = -divfaxr) in A r {x), 

(3.16) <\Vv\ p - 2 Vvn = on dA r (x) n (QU d N Q), 

v = on dA r (x) n 8dQ 

respectively. Notice in particular that equation (|3.15(l is well posed since its right hand side has 
zero mean value, because a is divergence free in f2, and a ■ n = on djsi£l so that 



div (ipa) dx = / div(ycr) dx = / a-ndTt 1 =0. 

n JA r JOA r 



We are now in a position to prove our minimality result. 



□ 



Proof of Theorem Vi.'A We claim that there exist at most m open balls B ri (x\), . . . , B rk (&&), k < 
m such that T C UiB ri (xi), 

n < CH 1 (B ri { Xi )n r), 

and 

B 2rj n B 2r;j (ay) = for all i ^ j, 

where C depends only on m. In fact let us consider the decomposition of T in its connected 
components, i.e., 

r := ri u • • ■ u r fc , 

with k < m. For all i = 1, . . . , k let B Si (jji) be an open ball with s, = ^(I^) and such that 

Tj Q B Si (yi). If the balls B 2si (yi) are disjoint, then the covering {B Si (2/i)}*=i ,fc satisfies the 

claim. Otherwise we proceed in this way. Let us consider 

fe 

S 1 :=|j5 2si (^), 

i=l 

and let fij, j = 1, . . . , fc < fc — 1 be its connected components. For all j, let _B Tj (zj) be an open 
ball with Tj — diam(Sj) such that 2?j C B Tj (zj). Again, if the balls B 2Tj (zj) are disjoint, then the 
covering {i? Tj g satisfies the claim. Otherwise we construct in a similar way as before 

the set B 2 which has at most k — 2 connected components. Clearly in at most m steps we come 
up with at most m balls satisfying the requirements of the claim. 

Since is Lipschitz, we have that if 7i 1 (r) is sufficiently small (depending on f2), we can assume 
that the balls {B ri (xj)}i = x 1 ... 1 fe intersecting dtt can be replaced by rectangles R Sj (yj) of the form 
(|3.1U|) satisfying l|3.11|l . centered at some point yj that falls into case (1), (2) or (3) of Lemma 13. 51 
More precisely, thare exist a constant C depending only on f2 and m, and there exist at most m 
open balls {B ri (xi)} and at most m rectangles {R Sj {yj)} defined in l|3.10|) with yj falling into case 
(1), (2) or (3) of Lemma 1331 such that B 2ri {xi) C Q, and 

r C Ui tj B ri (xi){JR Sj { yj ), 

(3.17) dtem{B ri (xi)) < CH^B^Xi) n T), diam(B Sj ( % )) < CH\R, •//.,: n T), 

B 2rj (a*) n B 2rj (ay) = 0, R 2St ( yi ) n i? 2Sj (%•) - 0, Bar, (a*) n R 2S] (yj) = for all 
Let y be a smooth function with 

< (p < 1, </? = on |J (B r< (xj) U fl^Gfe)) . V = 1 outside (J O^K) U R^M)) ■ 
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Let us denote with A ri (xi) and A Sj (yj) the sets Bi Ti {xi) \ B Ti (xi) and R2 S j(yj) \ RsAUj) respec- 
tively. Let 

77, S L q (A ri (x t );R 2 ) and r h E L q (A S] {y 3 );R 2 ) 

with q := p' = p/(p— 1) be the vector fields given by Lemmas 13.41 and 13.51 Let us consider 
■q £ L q {Q;R 2 ) defined as 

r\i in A n (xi), 
rjj in A s 
otherwise, 



7] := < 

and let us set 



r := ipa + rj, 

where a = df{x,Vu) is the stress of the elastic solution u. By construction we have that 



(3.18) t-\7v = for all v e A(T), 

Jn 

where A(T) is as in l|2.6[l . Moreover we have that for all i,j 



f \rj\ q dx < C f \\Ju\ p dx, 

J A h J A h 



where Ah denotes one of the A Ti (xi)'s or one of the A Sj (yj)'s, and C is a constant depending only 
on Q and the bulk energy density /. 

In view of H3.18(l we can put r in inequality H2.7|l getting 

(3.19) / [f(x, Vu) - f{x, Vu r )} dx < V / [r - a] ■ [df*(x, r) - df*(x, a)} dx. 

Jn j JA h 

Since /* is the convex conjugate of /, and since / satisfies the growth conditions (|2.3|l . we deduce 
that 

\t(x,0\ <c(ICI 9 + i) 

and 

\d c r( X ,o\<c(\cr i +i). 

We claim that for every Ah 



(3.20) 



/ [T-a}-[df*(x ) T)-df(x )( 7)}<Cdiaia(A h ) a , 



where C is independent of Ah, and depends only on m, and /. Here a > 1 is the exponent 
defining the weak singularities of u (see Definition ^. 1|) . From (|3.20|l and H3.19(l , taking into account 
H3.17(l and the fact that a > 1, we obtain that there exists I* depending only on fi, m, f and k 
such that for every L G JC m (fl) with H 1 (r) < I* we have 

/ [f{x,Vu) - f(x,Vu T )]dx < C Vdiam^)" < CH x (T) a < kH^T), 
J n h 
so that the minimality result holds. 

In order to conclude the proof, we need to show that claim (|3.20|) holds true. This can be seen 
making all the products and estimating each addend. Let us check the first one, the other ones 
being similar. We have 

r • df*{x, r)dx<C [ [\a\ + M) ■ (H^- 1 + + 1) dx. 



Then, in view of Lemmas 13.41 and 13.51 since u has uniformly weak singularities in fl we get for 
small enough 



/ \a\ q dx<C [ \\7u\ p dx < Cdiam(A ft ) a! , 

JA h JA h 

f \ri\ q dx<C f \\7u\ p dx <Cdmm(A h ) a , 

JA h JA h 
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\r}\\a\ q - 1 dx < (^J \n\ 9 <kcj q (^J \a\ q dx^" < Cdiam(A h ) a , 
\r]\ q - l \a\dx<(^j H q d x y f^J \a\ q d x y < Cdisan(A h ) a , 
[ \a\ dx < C ( [ \a\ q dx] " di&m{A h )i < Cdi&m(A h ) a , 

JA h \JA h J 



and 



/ \ri\dx<C[ [ \r/\ q dx) " diam{A h )p < Cdiam(^) Q . 
Summing up wc obtain that IpHH)) holds. □ 

Remark 3.6. Notice that the arguments of the previous proof also work in the case in which the 
elastic solution u has critical singularities, i.e. in the case a = 1, provided that there exists e > 
and 5 sufficiently small (depending only on O, m, f and k), such that 

Vit| p dx < Sr for every r < e. 

B r (x)nii 

Using the same arguments of the proof of Theorem 13.21 we can easily deduce the following 
localized version of the minimality result. 

Proposition 3.7. Let the elastic solution u of problem l|2.4|l have only uniformly weak singularities 
in A, where A is an open subset ofVL. Then there exists a critical length I* > depending only on 
A, m, f, and k such that for all T G K- m (fl) with T C A and 7i 1 (r) < I* we have 

[ f{x,X7u)dx< [ f[x,Vu T )dx + kH l {T), 

where ur is a minimum of l|2.4|l . 

Remark 3.8. (The case of strong singularities) Theorem 13.21 is false if the elastic solution u 
has strong singularities in fl, i.e., there exists x € such that 

(3.21) limsupi / \X7u\ p dx = +oo. 

r^o r JB r (x)nn 

In fact if such a point x exists, then the pair (v r , dB r (x) n f2), where v r is defined as 

u{y) in \ B r (x), 



(3.22) v r (y) 



o mnnB r ( x ) 



is energetically more convenient with respect to u for r small. Note also that this example needs 
the right hand-side of l|3.21|l to be just greater than 

The previous remark and Proposition ^ . 7l seem to suggest that in presence of strong singularities, 
energetically convenient small cracks prefer to stay near the singular points. Let us prove that 
this intuition is indeed true in the case in which u has only one point of strong singularity x. 

For every I > let Ti € /C m (f2) be such that (ur, , T;) minimizes the total energy (13.1(1 among all 
pairs («r,r) with T £ /C m (£!) and 7i 1 (r) < /. The existence of F; can be proved using the direct 
method of the Calculus of Variations in view of the lower semicontinuity of the T^-measurc with 
respect to Hausdorff converging sequences in K m (Vl) given by Gglab's Theorem (see for example 
OH for details). 

Notice that Yi ^ % for / small because the pair (vi,dBi), where vi is defined in 13.22fl . is 
energetically more convenient than the elastic solution. Moreover, for every r > 0, we have that 
F;nJ5,.(a;) ^ when I is small enough because otherwise, in view of Proposition l3~71 (up l , T;) would 
not be energetically more convenient than the elastic solution. So we deduce that the following 
proposition holds. 
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Proposition 3.9. Assume that x £ fi is a point of strong singularity for the elastic solution u, and 
that u has uniformly weak singularities inQ\ B r (x) for every r > 0. Then for every neighborhood 
U of x, if I is small enough we have Ti ^ 0, and TiHU =/= 0. 

Remark 3.10. (Singularities in materials) Theorem 13 . 21 and Remark 13.81 show that the quan- 
tity 

(3.23) f \Vu\ p dx, 

where u is a solution of 
(3.24) 



divd^f(x, Vu) = in fi 
u = ip on drjQ, 



determines the feasibility of the appearance of small cracks in ft when imposing a boundary datum 

Let us recall some known result concerning the behavior of (|3.23l) . As before, we assume that 
/ : £1 x R 2 — > K is a Carathcodory function satisfying 1|2.1|) and such that for a.e. x £ fl and for 
all £ £ K 2 

\\w<nx,z)<Mt\ p , 

where A, A > 0. 

By Theorem 13. 21 small cracks are energetically convenient if u has uniformly weak singularities 
in fl according to Definition 13.11 This is certainly the case if Vu has a summability sufficiently 
higher than p, as one can check by means of Holder inequality. In fact we have 

)|(«-p)/9 < / \Vu\ q dx) r 2 <- g - p V q . 
JB r (x) \JB r (x) J \Jn J 

We deduce that u has weak singularities in Q, if Vu £ L q (Vl\ K 2 ) with q > 2p. 

It is well known that if f2, ip and / are sufficiently regular, then the solution u of (|3.24(l is regular, 
so that u has only uniformly weak singularities in f2. However the assumption of regularity of 
/ with respect to the variable x is not suitable for applications to continuum mechanics, since 
discontinuity in the variable x models the important case of composite materials. 

Several papers in the literature address the issue of higher integrability properties of the gradient 
of solutions of l|3.24|l without continuity assumptions on x (see [T5] and references 

therein), namely properties of type Vu £ L* oc (^;]R 2 ) for some q > p. The restriction of the 
behavior of Vu to compactly contained open subsets of £1 is the natural price to pay in order to 
concentrate on properties depending only on the material (i.e., on the bulk energy density /) and 
not on the boundary datum. 

For example a regularity result by Caffarelli and Peral 7, Theorem C] implies that, under mild 
assumptions on /, if for all q > p there exists e — e{q) > such that if 

M p -H-9 i fM\<m p - 1 

(i.e., d^f is sufficiently near to the p-Laplacian operator), then u £ W^ q (fl). In particular if e is 
small enough, u has uniformly weak singularities on every compactly contained open subset of f2. 

In the linear case p = 2, the hypothesis of the regularity result by Caffarelli and Peral [J] reduces 
to the requirement that 

A/A < 1 + e 

for £ sufficiently small. However Leonetti and Nesi |15| found the optimal integrability exponent 
for Vu in terms of K :— (A/A) 1//2 , namely they proved that 

Vu G L s loc {Q) for every s < 2K/(K - 1). 

So in order to guarantee that u has uniformly weak singularities in compactly contained open 
subsets of f2, it suffices that A/ A < 4. 
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A famous example due to Meyers ^2] shows that if A/A = 4 strong singularities may appear 
inside the body. In fact, assuming that the origin belongs to SI we can consider 

A(x) := 1 /Kn ® n + Kt ® r, 

where n := A, r is obtained from n through a rotation of 90 degrees counterclockwise, and a ® b 
denotes the matrix with coefficients (a (g> b)ij = a%bj. Then it is easy to see that the solution of 

&iv(A(x)Vu) = in SI 
u(x) — xi on dSl 

is given by 

u(a;) := \x\ K ~ l xx. 

A simple computation shows that if K > 2, i.e., A/A > 4, the origin is a point of strong singularity 
for u. 

Remark 3.11. (The vectorial 2Z?-case) The minimality result given by Theorem 13 . 21 holds also 
in the case of R M -valued displacements provided that we choose them in the Deny-Lions space 

(3.25) L^ p (n\T;R M ) := {v e W^(n\T,R M ) : Vv £ L p (Sl \ T; M Mx2 )}. 

We have that W 1 *^ \ T,R M ) C L l > p (Q, \ T;K M ), and the two spaces are equal if SI \ T is 
sufficiently regular (for example an union of a finite number of Lipschitz domains). A notion of 
trace for functions in L 1,p (f2 \ T;R M ) near the points of doSl \ T is well defined, so that given 
ip £ W 1 ^^, M. AI ) and T £ K m (ST), the displacement Mr is a solution of the minimum problem 

(3.26) min {/ /(' T ' Vu ) dx : u e L hp {n\T;R M ),u = ip on <9 D f2\rj . 

(Notice that in this vectorial setting the maximum principle does not hold, so that only a control 
on the gradient is available and this is why Deny-Lions spaces are required). Since Definition 13. II 
relies only on the behavior of Vu, it turns out that the notion of uniformly weak singularities is 
well defined in the context of L 1:P -spaces. 

The minimality result in the vectorial setting follows because estimate Ij2.7ll still holds provided 
we set 

A(T) := {v e L^(n \ T; R M ) : v = on d D Q \ T}, 

and the constructions of Lemmas 13.41 and 13.51 can easily be adapted to the case of matrix valued 
vector fields. 

Remark 3.12. (The A-dimensional case) Let us consider the case f2 C with N > 3. The 
two dimensional setting is employed in the proof of Theorem 13 . 21 only to ensure the existence of a 
covering of the crack T which satisfies conditions H3.17J) . In the case T is connected, the covering 
condition can be reduced to the existence of an open set A (a ball if L is well inside SI, or a 
rectangle if it is near the boundary) such that T C A and 

diam(A) < CH^T), 

where C is a given constant. This inequality is not implied by connectedness in dimension N > 
3, because "needle-like" cracks can have very small 7i Ar_1 -measure and big diameter. So the 
machinery of the proof of Theorem 13.21 can be employed in dimension N > 3 if we restrict to a 
class of admissible cracks which excludes the elongated ones. Namely we can consider the family 

(3.27) /C C (H) := {r C U : T is closed and diam(r) < CTi w_1 (r)}, 

where C is a given constant. The displacement Mr relative to a crack T S K, (SI) and the boundary 
datum tp £ W 1,p (Sl,M. ), with M > 1, is given again by problem 

(3.28) min !j f(x, Vm) dx : u £ L x ' p (n \ T; R M ), u = ip on d D Sl \ r| , 
where the space L 1 *^ \ T;R M ) is defined in (|3~25j) . 
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The notion of uniformly weak singularities can be rephrased in the TV-dimensional setting in 
the following way: we say that u € L 1,p (f2; R^) has uniformly weak singularities in Vl if there exist 
constants N — 1 < a < N and C > such that for every x € f2 and for r small enough 

(3.29) [ \X7u\ p dx <Cr a . 

JB r {x)nn 

The intuitive meaning of condition (|3.29l) is the same as in the two-dimensional setting, namely 
that the bulk energy of u inside a ball of center x and radius r is asymptotically negligible for 
r — ► and uniformly in x with respect to the surface of the ball. 

4. The minimality result in linearized elasticity 

We briefly show in this section that the results obtained up to now are also true if the bulk 
energy f(x,\7u) is replaced with a quadratic linearized elasticity energy 

C(x)e{u) : e(u) = C it j fkf i(x)e{u) it je{u) k ,h 

that is a positive-definite quadratic form of the symmetrized gradient 

e(u) := (Vu + (Vu) T )/2 

of the displacement u : fl — > M. N . As before, the main conclusions will be drawn in dimension 
2 (Theorem 13.21 or rather, since u is vectorial, the result of Remark 13 . 1 1|> . whereas in higher 
dimension the same restrictions of Remark 13 . 1 21 will apply. 

Let us consider in fl a measurable function C, such that for any x € fi, C (x) is a N x N x N x N 
4th-order tensor that defines a positive-definite quadratic form on the vector space of symmetric 
N x N matrices, that we denote by S NxN . We assume that for any £ € gNxN anc j a e x g ^ 
holds 

(4-1) X\e < C(a:)£:£ < A|£| 2 , 

where £ : 77 = Tr(£i] T ) = £i,jT]ij and |£| 2 = £ : £ is the standard Euclidean (Frobenius) norm. 

Given T C f2 a compact one-dimensional fracture, the space of admissible displacements with 
finite energy will be the space of measurable displacements u : f2 — > whose symmetrized 
distributional gradient in CI \ T, denoted by e(u), is in L 2 (Q, \ T; S NxN ), and that satisfy in some 
sense u = ip on Thanks to Korns inequality , it is known that such a displacement belongs 

in fact to £/jQ C (fi \ T), and since we have assumed that the boundary of £1 is Lipschitz, we also 
have that u E 1 (il n B) for any ball B with B n T = 0, so that the trace of it on f2 \ T is well 
defined. As in Section 13. Ill we introduce the space 

LD^{Cl\T) :={ueHUn\T;TSL N ) : e(u) G L 2 (fi \ T;S NxN )} 

and for T a closed subset of 17, the displacement ur is given by 

min j^" C(x)e(u) : e(u) dx : u £ LD h2 (Q\T),u = ip on a D fi\r| . 

We denote by u the solution in the non-cracked domain. The set of admissible variations is now 

^(T) := {v e LD h2 (fl \ T) : v = on d D d \ V}. 
The proof of Theorem 12.11 can be reproduced in this situation, yielding the estimate (we assume 

|r| = 0) 

f [C(x)e(u) : e(u) - C{x)e(u r ) : e(« r )] dx < 2 f [t - a] : [C(:t) _1 (t - a)] dx, 
Jn Jn 

with <j{x) := C(x)e{u){x) and r € i 2 (^2; S NxN ) is any stress field compatible with the variations 
in A(T), that is, such that 

(4.2) / r : e(v) dx = 

Jn\r 
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for any v £ A(T). The last estimate, together with Q4.1|) . yields 

(4.3) / [C(x)e(u) : e(u) - C{x)e(u r ) : e(u r )] dx < \ I |r - a\ 2 dx , 
Jn * Jn 

for any r satisfying (|4.2H . 

In order to prove Theorem 13.21 in this setting, we do exactly the same construction as before. 
We build a r from a by letting r = a in f2 except in a finite union of balls or rectangles Bi Ti {x,i) 
or i?2s (%•) (cf the proof of Theorem 13 . 21 dp . [TU1 and following). Inside the smaller balls/rectangles 
B Ti (xi) and R Sj (yj), we choose r = 0, and in each crown A rj (a)i) = B2 ri {xi)\B ri (xi) or v4 Sj (j/j) = 
Ifasi illj) \ R Sj {yj), t is of the form: r = <^<t + r/ where (/? is the appropriate cut-off function. 

Again, to achieve l|4.2|l . one needs to choose r\ in an appropriate way. An additional difficulty 
here follows from the fact that rj(x) has to be almost everywhere a N x N symmetric matrix. 
In order to find a suitable r/, one needs to replace problems (|3.7|) . (|3.14|) . 13.15|) or l|3.16|l by the 
appropriate variant. 

The right way to do is obviously to solve the vectorial equation 

(4.4) dive(w) = — div ((per) 

in the appropriate domain, and to replace the Neumann boundary condition, when present, with 
the corresponding condition e(v) ■ n = 0. We then set r; = e(v) in each crown. 

We get an estimate on J n |r — tr| 2 dx from standard estimates on e(v), that will follow from 
appropriate Poincare-Korn inequalities. For instance, one shows (see Appendix [B] for details) that 

(4.5) / \v\ 2 dx < Cr 2 [ \e{v)\ 2 dx 

J A r (x) J A r (x) 

for any v € if 1 (A- (a); l") with 



v(y) dy = and / y x v(y) dy = , or v = on <9of2 n dA r (x) . 

A r (x) J A r (x) J 

The first set of conditions ensure that the "average rigid motion" of v vanishes in A r (x). Rigid 
motions, of the form x i— ► Kx+p with K antisymmetric, arc the kernel of the symmetrized gradient 
(in any connected domain). 

Multiplying l|4.4ll by v and integrating by parts, we find that 

\e(v)\ 2 dx = - I a : (Vip ® v(x)) dx 

A r {x) J A r (x) 

< \W\\L 2 (A r (x))\\^f\\L^(A r (x)) (yCr\\e(v)\\ L 2 {Ar{x)) ^J , 
where we have used (|4.5|) . Since || Vt^H^oo^^ ^ < 2/r and rj = e(v), we deduce that 

|77(ar)| 2 ^ < C I \cr(x)\ 2 dx 

I A r (x) J A r (x) 

for some constant C that does not depend on x or r. We conclude as in the proof of Theorem 13. 21 

5. Qualitative properties of crack initiation 

In this Section we use the results of Section|21to address the problem of crack initiation in elastic 
bodies. We restrict our analysis to the case of antiplane elasticity. In view of the minimality result 
of Section the same conclusions hold also for the case of planar linearized elasticity. 

First of all we consider the classical Griffith's theory of quasistatic crack propagation, and we 
prove that it cannot explain the formation of a crack in an elastic body Q without singularities 
within the class /C m (f2) of closed sets with a finite number of connected components, which is 
much richer than the family of smooth curves usually considered in the literature. 

In the second part of the section, we remove the assumption, implicit in Griffith's theory, that 
the crack growth is progressive, namely that the length of the crack is continuous in time. Inspired 
by the variational theory of quasistatic crack evolution proposed by Francfort and Marigo in |13| . 
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we replace the classical Griffith's equilibrium condition with a static equilibrium condition and an 
energy balance. The static equilibrium condition is a unilateral minimality property which states 
that, during the crack evolution, the total energy is minimal among all configurations with larger 
cracks (so that discontinuities of the crack's length are allowed). The energy balance requires that 
the total energy of the system evolves in relation with the power of external loads in such a way 
that no dissipation occurs. 

Within this framework, we prove that a crack appears immediately at a point of strong singu- 
larity for the body. Moreover we prove that if the body has uniformly weak singularities, then it 
deforms elastically until a critical time ti after which a "big" crack T(t) appears. These results 
have been established by Francfort and Marigo in Proposition 4.19, point (ii)] under the as- 
sumptions that the crack T(t) is union of m fixed curves {7i(t)}i=i,..., m which can be parametrized 
by arc length. Thank to our local minimality result ( Theorem 13. 2f> . we prove these facts removing 
the restrictions on the path of the crack. 

The mathematical setting we consider is that of Section [21 Namely £1 is a bounded Lipschitz 
open set in R 2 , do®* Q dVl is open in the relative topology, and d^Cl := d£l \ do^ is composed 
of a finite number of connected components. The family of admissible cracks is given by the class 
/C m (f2) of closed subsets of with at most m connected components and with finite length, while 
the class of admissible displacements relative to a crack T is given by W 1,p (Vl\T) with p +oo[. 
The total energy is given by 

(5.1) £(u,T):= / f(x,Vu)dx + kH l {T), 

where / is a Caratheodory function satisfying (|2.1[l . (|2.2I) and (|2.3|) . and k > 0. 



5.1. Crack initiation and Griffith's theory. Let To be a crack inside O of length lo, and 
suppose that a boundary displacement ifi is assigned on 8dQ \ T. According to Griffith theory, Tq 
is in equilibrium if, taken any family of increasing cracks Ti containing To with length lo + l, then 

/, s ,. W(l )-W(l + l) 
(5.2) hmsup < k, 

where W(Iq) and W(lo + 1) denote the bulk energy of the displacements ur and wr ; associated to 
the boundary datum ip and the cracks To and Ti respectively, and k represents the toughness of 
the material. Moreover, during a quasistatic crack evolution, if To propagates along the T/, then 
(|5.2|) holds with equality. 

Let us prove that the rate of energy release that appear in the left hand side of l|5.2|l is zero in 
the case in which F = and the elastic solution u relative to the boundary displacement tp has 
uniformly weak singularities. This means that the elastic configuration is always in equilibrium 
according to Griffith's theory, and moreover that a quasistatic crack evolution which begins in the 
elastic configuration remains at all subsequent times in the elastic regime, i.e., Griffith's theory 
cannot explain crack initiation. 

To this aim for every I > 0, let us set 

W(Z) := inf jy f(x,Vu T )dx : T <E K m (fl), U 1 (F) < zj , 

where u-p £ W 1,p (tt \ T) denotes the displacement associated to F and the boundary datum ip. 
Notice that we have clearly W(0) = J n f(x, Vu) dx, where u is the elastic displacement relative to 

The following proposition holds. 
Proposition 5.1. Let us assume that the hypothesis of Theorem \H .°A are fulfilled. Then we have 

hm MO) -MQ = 
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Proof. For every k > 0, by Theorem 13.21 we have that for I small enough and T G /C m (f2) with 



/ f{x,\7u)dx< [ /(x,Vu r ) + fcH 1 (r). 



We deduce that 

w(o)-w(Q - 

hmsup : < k. 

Since fc is arbitrary, and since W(l) < W(0), we conclude that the result holds. □ 

Remark 5.2. (The case of strong singularities) If the elastic solution u has a strong singu- 
larity at x € f2, then by Remark 13.81 we have that 

,. W(0)-WQ) 

hm — — — = +oo, 

z^o+ I 

so that the elastic configuration is not in equilibrium in the framework of Griffith's theory. 

5.2. Crack initiation in variational theories of crack propagation. As explained at the 
beginning of the section, we consider now irreversible quasistatic crack evolutions governed by a 
static equilibrium condition and an energy balance. More precisely if ip(t) is a time dependent 
boundary displacement, and u(t),T(t) are the displacement and the crack at time t relative to 
ip(t), we assume that the pair (u(t),T(t)) satisfies the following properties: 

(a) Irreversibility: T(t) is increasing in time, i.e., T(tx) C T(i 2 ) for all < t\ < t 2 < T; 

(b) Static equilibrium: if t > 0, £(u(t),T(t)) < £(u,H) for all cracks H such that U s <tr(s) C 
H and all displacements v : \ H — > R with v = ip(t) on <9of2 \ H; 

(c) Energy balance: the total energy £(u(t),T(t)) is absolutely continuous in time, and it 
satisfies 

£(u(t),T(t)) = £(u(0),r(0))+ / / df(x,Vu(T))Vip(T)dxdT. 

Jo Jn 

Condition (a) stands for the irreversibility of the evolution: the crack can only increase in time, 
i.e., no healing processes are admitted. Condition (b) asserts that the pair (u(t), T(t)) is a unilateral 
minimizer of the total energy, i.e., it is a minimum among all configuration with larger cracks. 
In particular u(t) is the elastic deformation relative to the boundary datum ip(t) in the domain 
f2 \ r(i), i.e., u(t) satisfies equation l|2.5l) with T — T(t) and ij; = 4>(t). Finally, notice that under 
suitable regularity assumptions on u{t) and T(t), condition (c) can be rewritten as 



f(x,Vu(T)) ■ ijj^dH 1 ^) dr. 



(5.3) / f(x,Vu(t))dx - / /(x^uis^dx + k^iK^XKis)) 
Jn Jn 

d_ 

d D n\K(T) dn 

Therefore the energy balance condition states that the sum of the variation of the bulk energy and 
of the dissipation due to the creation of a new crack is equal to the work inserted in the system 
by the boundary datum ip. We refer the reader to the paper by Francfort and Marigo 13 for 
further details on crack evolutions satisfying conditions (a), (b) and (c) (see also Mielke ^H] for a 
connection with the theory of rate- independent processes). 

In order to treat the problem of crack initiation, we consider as in J3j the case in which / is 
p-homogeneous in the gradient, i.e., for all x e £ e K 2 and t > 

f(x,tO=t p f(x,0- 

We consider a time dependent boundary displacement of the form t — > tip, where t G [0, T] and 
ip <G T / F 1,p (il) n L°°(f2) is a given function. We refer the reader to the paper by Dal Maso and 
Toader |U] for the existence of a quasistatic crack evolution satisfying (a), (6) and (c) within this 
mathematical setting. 

Since we are dealing with a crack initiation problem, and since "0(0) = 0, we assume that 
(w(0),r(0)) = (0,0). Notice that if?; denotes the elastic displacement associated to the boundary 
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datum ip, then tv is the elastic displacement associated to tip. Then from the static equilibrium 
condition, comparing (u(i),T(t)) with (tv,T(t)), we have that for all t g]0, T] 



(5.4) / f(x,Vu{t))dx <t p / f(x,Vv)dx. 

Jn Jn 

Finally, since we can replace tip by tu in the energy balance condition, we can write 

(5.5) / f{x,Vu{t))dx + kH 1 {T{t)) = / / df{x,Vu{r)) ■ Vvdxdr. 
Jn Jo Jn 

This implies that 

(5.6) \ f(x,Vu(t))dx + kH l (T(t))<t p [ f{x,X7v)dx. 
Jn Jn 

In fact, from the inequality 

xy<f(x) + r(y), 

and taking into account the p-homogeneity of /, we obtain the following Holder type inequality 

(5.7) [ zw< p^q 1 ^ ( [ f(x, z) dx) ^ ( [ f*(x, w) dx) ^ , 
Jn \Jn ) \Jn J 

where q — p' := p/(p — 1). Taking z = Vw and w = df(x, Vw(t)) we get 

i i 

f df{x,Vu{T))Vvdx<p 1 / p q 1 ' q [ f f{x,Vv)dx ) ( f f*{df(x,Vu{T)))dx) q 
Jn \Jn J \Jn J 

i_ i 

= P 1/p q 1/q (J n f(x, V«) dxj " (p - 1)1 U f(x, Vu(r)) dx) ' . 

In view of (|5.4|l we have 

/ df(x, Vu(r))Vw dx < T*p 1/p q 1/q (p- l)i / f{x,S7v)dx. 
Jn Jn 

Integrating from to t we obtain 

df(x, X7u(T))X7vdx <p xlv q 1 l q p- 1 {p- l)*t p ( f{x,Vv)dx. 



o Jn 



Since 

pVp g V9p-i(p_i)i = 

by (|5.5|l we conclude that (|5.6|l holds. 

Notice that we can rescale l|5.6|l obtaining for t small 

(5.8) / f{x,\Jv{t))dx + H l {T{t)) < ( f(x,Vv)dx 



where v(t) := \u(t) is the displacement associated to T(t) and ip. 

As noticed by Francfort and Marigo in |13j . if T is large enough, a crack will appear during the 
evolution, i.e., T(s) ^ for some s G]0, T[. In fact if T is such that 

(5.9) H 1 (d D n)<T p I f(x,Vv)dx, 

Jn 

then we get that creating a crack along 8d^ is more convenient that deforming elastically. We 
are now in a position to state the first crack initiation result. 

Theorem 5.3. Let us assume thatT satisfies l|5.9|l . and let us suppose that the elastic displacement 
v associated to the boundary datum ip has uniformly weak singularities in Q,, i.e., it satisfies l|3.3|l . 
Then the crack initiation is brutal, i.e., there exists a positive time U £]0,T] such that T(t) = 
for every t < U, and 7i 1 (r(t)) > I* for all t G]U, T] for some I* > depending only on fl, f , k, m 
and ip. 
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Proof. Notice that by f5U|l we have H x (r(t)) -> as i -> 0. By Theorem we get that for t 
small the elastic solution v is energetically convenient with respect to (v(t),T(t)), where v(t) is 
the displacement associated to T(t) and ip. But this is against i|5.8|) unless 7i 1 (r(t)) = 0. In view 
of (EH we deduce that there exists U g]0, T[ such that T(t) = for every t < U and W 1 (r(t)) > 
for allf e]tf,T]. 

In order to prove that W 1 (r(i)) > I* for all t e]ti,T] and for some /* > 0, notice that by (|5.6I) 
we deduce that 

f /(x,Vu(t))rfa; + fcT- p W 1 (r(t)) < / f(x,Vv)dx. 
Jn Jn 

Then by Theorem 13. 21 we deduce that 

lim mfT^CTft)) > I* 

for some positive constant I* depending only on f2, /, kT~ p , m and ip. The proof is thus concluded. 

□ 

Let us now study the crack initiation in the case in which the elastic displacement u associated 
to the boundary datum ip has strong singularities. We recall that a point x € £1 is a point of 
strong singularity for v if 

limsup- / | Vw| p dx = +oo. 

r^o r JB r (x)nn 

It is well expected that a during a loading process, a crack will appear at a point of strong 
singularity. The following theorem establishes this fact for a general quasistatic crack evolution 
satisfying properties (a), (b) and (c). 

Theorem 5.4. Let us suppose that the elastic displacement v associated to the boundary datum 
ip has a strong singularity at x £ and that v has uniformly weak singularities in O \ B r (x) for 
every r > 0. Then we have that 7i 1 (r(t)) > for all t G]0, T], and the crack starts at the point 
{x}, i.e., 

xe f|r(t). 

t>0 

Moreover the crack departs with zero speed, i.e., 

(5.10) hm^™1^0. 

Proof. In view of Remark 13. 81 and since x is a point of strong singularity for v, we have that for a 
positive time t > the elastic displacement to relative to the boundary datum tip cannot satisfy 
condition (b). As a consequence, we deduce that 7i 1 (r(t)) > for all t s]0,T]. 
Let us come to the properties of T at time t = 0. By (|5.6|l we have 

limft x (r(t)) = 0, 

so that 

f|r(t) = {»i,... ,»h}, 

t>0 

with /i < m. Let us suppose by contradiction that yi ^ x for all i — l,...,h. Then there 
exists r > such that for t small enough we have T(t) C f2 \ B r (x). Since u has uniformly weak 
singularities in A := O \ B r (x), and in view of Proposition 18.71 we have that inequality Ij5.8|l 
implies that 7i 1 (r(i)) = for t small, which is a contradiction. 
Finally, in order to prove (|5.10() . we rescale l|5.6|l obtaining 

/ fix^vit^dx + ^-H^Tit)) < f f{x,Vv)dx, 
Jn 1 Jn 

where v(t) is the displacement associated to T(t) and ip. We deduce that 7i 1 (r(i)) < Ct p for some 
constant C, and so (|5.1U|) easily follows. □ 
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Appendix A. The two dimensional SBV case 

The aim of this appendix is to prove a minimality result in the lines of Theorem 13.21 which 
does not require an a priori bound on the number of the connected components of the admissible 
cracks. Small cracks are still not energetically convenient if the gradient of the elastic solution of 
problem Ij2.4|l and the related stress are continuous in f2. This condition excludes however that 
the elastic configuration presents (weak) singularities. 

In order to make the mathematical setting of this section precise, we need to recall some facts 
about rectifiable sets and the functional space SBV of special functions with bounded variation. 
We refer the reader to 0] for a complete treatment of these subjects. 

A set r C M. N is rectifiable if there exists N C T with H N ~ x (N Q ) — 0, and a sequence (Mj) ig ^ 
of C 1 -submanifolds of M. N such that 

T\N c[jM t . 

zgN 

For every iem Mi , we define the normal to T at x as Um^ (x) . It turns out that the normal is 
well defined (up to the sign) for 7i N ~ l -&.e. x E T. 

Let U C M. N be an open bounded set with Lipschitz boundary. SBV(U) is the set of functions 
u G ^(U) such that the distributional derivative Du is a Radon measure which, for every open 
set AC. U, can be represented as 

Du(A) = I Vudx+ I [u]{x)vdH N ~ l {x), 

J A JAnS(u) 

where Vu is the approximate differential of u, S(u) is the set of jump of u (which is a rectifiable 
set), v(x) is the normal to S(u) at x, and is the jump of u at x. 

For every +oo[ we set 

SBV P (U) := {u G SBV(U) : Vu G L P (U,R N ), H N ~ 1 (S(u)) < +oo}. 

If u G SBV(U), then u admits a trace j(u) on dU which is characterized by the relation (see 
Theorem 3.87]) 

lim r~ N I \u(y) - -y(u)(x)\ dy = for H N ~ x -a,.e. x G dU. 

JnnB r (x) 

We will denote the trace 7(1*) on dU again by u. If T C U is rectifiable and oriented by a normal 
vector field n, then we can define the traces 7^ (u) and 7p (u) of u G SBV{U) on T (see 0] 
Theorem 3.77]) which are characterized by the relations 

lim r~ N I \u(y) - ^(u)(x)\ dy = for H N ~ l -a.e. x G T, 

where Bf{x) := {y G B r (x) : (y — x) ■ n ^ 0}. 

A set ECU has finite perimeter in U if the characteristic function 1# belongs to SBV(U). 
We denote by d*E the set of jumps of 

Let us now come to our problem of local minimality. Let fl C M 2 be open, connected and 
with Lipschitz boundary. Let Soil C dQ be open in the relative topology, and let 9/vf2 := 
dft\dfl. Let f(x,£) be a Caratheodory function satisfying l|2.1|l . (|2.2(l and (j2.3(l . Let us moreover 
assume that the boundary displacement on <9dJ7 is given by the trace of a continuous function 
i/j G C (D) n W l ' p (Q). We denote by u the elastic solution relative to ip, namely the solution to 
the problem 

min I J f(x, Vu) dx : u G W rl ' p (fi), u = ^ on d D n^ . 

The class of admissible cracks we consider is 

(A.l) TZ(U) := {T C U : F is 1-rectifiable and H^T) < +00}. 

Let us come to the class of admissible displacements relative to a crack T and to the boundary 
displacement ip. Since T is not supposed to be closed, the Sobolev space W 1 ' P {$1 \ T) is not well 
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defined. So we consider as class of admissible displacements the functions u G SBV p (Vl) such that 
S^{u) C r, where 

S*{u) := S{u) (J{x£ d D Q : u{x) 7^ i/j(x)}, 

and the inequality on <9d57 is intended for the traces. Notice that if T is closed then u S W ,p (il\T), 
and u — ip on <9dJ1 \ I\ 

The displacement ur € 5.BV p (f2) PI L°°(f2) associated to T and ip is a solution of the minimum 
problem 

(A.2) min | J f(x, Vit) dx : u € SW^Q), S^{u) C T 

The proof of the existence of ur is standard: it relies on Ambrosio's compactness and lower 
semicontinuity theorem [2], together with a truncation argument. 
The main result of the section is the following one. 

Theorem A.l. Let u be the elastic displacement relative to £ W ,p (fi) n C (0), and fei us 
assume that u satisfies 

(A.3) V«eC°(0;K 2 ) and a := d^f(x,Vu) G C°(n;M 2 ). 

T/ien i/iere exists a critical length I* > depending only on fl, f and ip such that for all T G 1Z(CI) 
with Ti. 1 ^) < I* we have 

[ f{x,\7u)dx< I f(x ) Vu r )dx + kH 1 {r). 
Jn Jn 

In order to prove Proposition lA.il we need the following lemma. 
Lemma A.2. For every T G 7£(f2) we have 

(A.4) [ [f(x,Vu) - f(x,Vu r )]dx< [ a ■ n{u+ - Up) dH 1 , 

Jn Jr 

where a is the stress of the elastic displacement u defined in (|A.3|) . uy is a minimum of (|A.2f) . 
and Mp are £/ie traces of ur on T. 

Proof. By the convexity of / we have 

(A.5) / [f(x, Vu) - f(x, Vu r )] dx < / df(x, V«)(Vw - Vu r ) dx = / ct(Vu - Vw r ) cte. 

it! it! JO 

We can assume that ip is defined on R 2 , i.e., ip G W 1 ' P (R 2 ) n C°(R 2 ). Let 5 be a ball centred at 
such that £1 C B. Let us set Q' := B\ c^O. We can extend u and ur to Q! setting u = ur = ip 
on B \ f2. Let us consider u n G C 1 (ri / ) with it„ = ip on B \ f2 and such that 

Du n — > Du strictly in the sense of measures, 

that is (see H Theorem 3.9]) 

lim / <pdDu n = I <pdDu forall( y 3GC (fi')ni oo (f7'). 

n^ + oo J n , J n , 

Since u is a minimum for problem (| A. 2|> with T = 0, and u n = ip on c?^)Sl, then we have that 
u — u n is an admissible variation for it so that for all n G N 

(A.6) 0= / a(Vu - Vu„) dx = / a dD(u~ u n ). 

Jn Jn 

Since by assumption a G C°(fl), we can extend cr to f£' in such a way that a G C°(f2') n L°°(fl'). 
Notice that since u = u„ = ip on B \ Q., from ifO)) we have 

a dD(u n — u) = 0. 
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Then by strict convergence we deduce 

(A. 7) = lim / adD(u-u n )= / aVudx- / a dD 
™^°° Ja.> Jn< Jn> 



/ it(Vii — V«r) dx — / g ■ n(u^ — u r ) dTL 1 . 
JiV Js(ur) 



'S{u r ) 

Since S(ur) C r, ujt = itp on T \ 5(itr), and Vu = Vur = VV> on -B \ fi, by (|A.7(I we deduce 

/ cr • n(iip — Mp ) dTL 1 — / cr(Vu — Vur) da; = / <r(X7u — X7ur) dx, 
Jr Jfi' Jn 

so that, in view of (TOl . we have that (|A.4(I follows. D 

We are now in a position to prove the minimality result in the SBV context. 
Proof of Provosition \A.lj Let us consider I > 0, and let T £ TZ(fl) be a minimum for the functional 

jP(r) := / fix^urjdx + kH 1 ^) 



among the cracks L E 1Z(H) such that 7i 1 (r) < I. The existence of such a T follows minimizing 
the functional 

F(u) := / /(x,Vu r )^ + fcW 1 (S" /i (u)) 



among all u S SBV(Cl) with ^(^(u)) < /, and choosing P := S'(u). Proposition lA~Tl will be 
proved if we show that r = for I small enough. 

Let us establish some qualitative properties of T that are useful in the proof. Let E C S7 be a 
set with finite perimeter such that d*E C r. We can assume that 

(A.8) u r = on E. 

In fact we have that 

{ur in f2 \ E 
in E 

belongs to SBV(il), and it is an admissible displacement for T and ip (see 0J Theorem 3.84]) with 

f(x,Wur)dx< / f(x,Vur)dx. 



We conclude that fir is a minimum energy displacement relative to T and ip, so that l|A.8J) follows. 

Let i?i and E2 be two sets in fi with finite perimeter such that d*E\ C T and d*E2 C T. We 
have that 

(A.9) H 1 (^in9*E 2 )=0. 

In fact, if by contradiction TL 1 (d*Ei n d*E2) > 0, let us consider the crack 

f :=r\(a , B 1 na , B 2 ). 

Notice that S(u r ) Q f. In fact S(u r ) Q L, and u r = on Ex U _E 2 in view of (|A. 8|) so that 

(5(« r ) n n <9*£ 2 )) = 0. 

Since 7i 1 (L) < 7i 1 (L) < £ and Mr is an admissible displacement for V, we get that !F(T) < .F(r), 
and this is against the minimality of T. Then i|A.9|) holds. 

As a consequence of (|A.9|) . there exists a family (E^g^ of sets with finite perimeter in fi such 
that d*Ei C T, H 1 (d*E h n d*E k ) = for every h ^ k, and such that if E is a set with finite 
perimeter in fi such that 9*E C L, then E = Eh for some /i € N. Let us set 

00 

(A.10) T* := |J 8*E l7 

i=0 
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and let us assume that n denotes the outward normal to E^ . Since the stress a is a divergence free 
vector field, in view of the Generalized Gauss-Green formula for sets with finite perimeter (see 0] 
Theorem 3.36]) we have for all i G N 

(A.ll) / o-ndH l = / divadx = 0. 

In order to prove our minimality result, by contradiction let us assume that there exists lh — > 
and such that setting Uh :— ur h we have 

(A.12) / f(x,Vu h )dx + kH 1 (T h ) < [ f{x,Vu)dx. 

Jn Jn 

By Ambrosio's lower semicontinuity Theorem [2] and by l)A.12|) . we deduce that for every open 
set A C Q 

(A.13) lim / f{x,Vu h )dx = \ f{x,Vu)dx. 

h-^+oo J A J A 

Let (E^)i e fi and be the sets associated to Th described above. We claim that 
(A. 14) \[ u h]\ '■= \u£ — u7 | — * uniformly on Th, \ T^ as h — > +oo 

and 

(A. 15) ess-sup a * B h ( u h) ~ ess ~' m ^d*E h ( u h) ~ > ^ uniformly in i as h — > +oo. 

In view of (|A.14() and (|A.15() . the proof of the proposition is readily concluded. In fact, given 
e > 0, and choosing h so large that 

IM<£ on r h \n 

and for all i e N 

ess-sup a . B h(u£) — ess-'mi d , E h(u^) < e, 
by Lemma IA.2I and in view also of i|A.8() and of IjA.lip we have that 

/ [f{x, Vu) - f(x, Vu h )] dx< [ a- n{u+ - u^) dU l 
Jn Jr h 

<[ cr-n(u+ -u^)dH l + V / a-nu+dH 1 

Jr h \r> h l= o J d'E^ 

oo „ oo 

< eMooH 1 (T h \ T* h ) + V(ess-inf a , £flU +) / a ■ ndU 1 + e|M|oo V H l {d*E%) 

= eiHloo (n 1 (r fc \ r%) + n 1 (r;j) = eM^n 1 (r h ) , 

and this is against l|A.12|) . 

In order to conclude the proof, we have to prove the claims l|A.14|) and (|A.15|I . Let us consider 
(|A.14|I . the proof of the other claim being similar. Let us assume that there exists 8 > and 
Xh € Th \ T^ with 

(A.16) \[u h ](x h )\>5>0. 

Up to a subsequence we have Xh — » x £ f2. Let us assume that x E fl. For h large enough, and for 
r small we have B r (xh) C fl. Notice that for every r such that dB r (xh) H Th = 0, we have that 
Uh G W 1 ' p (dB r (xh))- Moreover by the maximum principle we have that 

(A. 17) max Uh — min Uh = Mh — mh > 5. 

dB r {x h ) dB r (x h ) 

In fact otherwise, we can consider Uh defined as 



u h 



Uh outside B r (xh), 

m&x{mm{uh, Mh},rrih} inside B r [x h)- 
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Since x^ G T;, \ T^, and in view of (|A.16(1 we deduce that 

\{u h ?u h }\\jE?\>a 

so that 



/ f(x,Vu h )dx< \ f(x,Vuh)dx 
Jn Jn 



which is against the minimality of Uh- Then (|A.17|I holds. 

Let dftUh denote the angular derivative of Uh, i.e., d$Uh :— S}'Uh( x h + r cos??, yh + rsin$). 
Setting C h := {s G [0, r] : H 1 {dB s (x h ) n Ift) = 0}, by (|XT7)l we have for every r e C h 



2tt 

\d$u h \ p dd > 2ir8 p . 



Since \Ch\ > r — H}(Th), we deduce that 



(A.18) / \Vu h \ p dx> / / s 1 - p \d i} u h \ p d'dds 

J B r {x h ) JO JO 

■ / s 1 - p \d^u h \ p ddds>2Tr8 p s 1 - p ds>27r6 p s^~ p ds. 

Jc h Jo Jc h J-HHTh) 



Let us distinguish two cases, namely p > 2 and 1 < p < 2. If p > 2, choosing r = 27i 1 (r/ l ) we 
obtain 

lim inf / f(x, Viifc) dx > 2%S P In 2. 

h^+oo J Brh(xh) 

But this is against l)A.13|l : in fact for all r such that B r (x) C ft by (1 A. 13(1 we have 

/ f(x,Vu)dx= lim / f(x,Vu h )dx> lim / /(x, Vu h ) dx > 2nS p In 2, 

ifl r (s) h ^+°°JB T {x) h ^+°° J B Th {x h ) 

and this gives a contradiction for r small enough. If 1 < p < 2, then we have 

lim inf / f{x,Vu h )dx > lim inf-— (r 2 - p - H 1 (T h ) 2 - p ) = - 

h^+oc J Br(xh) h^ + co2-p y 2-p 

from which by (|A.13J| we deduce that 

(A.19) f f(x,Vu)dx>Cr 2 - p 

J B r (x) 

for some C > 0. Since Vu G C°(fl;R 2 ) and / satisfies (j2.3j) we get that 

/ f{x, S7u)dx < Cr 2 

JbAx) 



r 2~P 



for some C > 0, which together with (|A.19|) gives a contradiction. 

The case in which x € dQ can be treated almost in the same way as the case x G f2. In fact it 
is sufficient to choose r so small that 

max ib — min ib < 5, 

B r (x h )nd D n B r {x h )nd D Sl 

and to take into account the fact that there exists a constant C depending only on f2 such that 

H l (dB rh (x h )nfl) > Cr h . 

In this way, integrations involved in l(A.18|l can be performed on a set of angle d which has 
a positive measure uniformly bounded from below, and the contradiction follows by the same 
arguments used above. □ 
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Appendix B. Uniform Poincare and Poincare-Korn inequalities 

In this section, we show that very basic arguments lead to the uniform Poincare and Poincare- 
Korn inequalities that are needed respectively in Sections [3] and 0] 

To simplify, we will only consider the case of iV-dimcnsional domains of the form 

Q f = {x = (x',x N ) ER N : < Xi < l,i = 1,.. .,N- 1,0 < x N < f(x')} 

where / : Q' —> [1, M] is a L-Lipschitz function. Here Q' is the (N— l)-dimensional cube (0, l)^ -1 
and L > and M > 1 are fixed constants. The adaption of the argument that we will present here 
to the "real" cases that are useful in the paper is straightforward. Also, for simplicity, we consider 
here the "linear" case p = 2. However, the proofs would be identical with any other exponent 
p£ (l,+oo). 

With a slight abuse in the notation we also identify Q' with the base of Qf, that is, the subset 
(0, l)^- 1 x {0} of dQ f . We show that the following result holds: 

Proposition B.l. There exists a constant C > depending only on L and M such that 

(i) For any u £ H 1 (Qf) with u = on Q' , ||m||l 2 (q / ) < C||^ u IIl 2 (Q/)/ 

(ii) For any u £ H 1 (Q } ) with / u{x) dx = ; ||w|| L 2( Q/) < C|| Vu|| i 2 (Q/ ); 



(iii) For any u £ H^Qf,^) with u = 0on Q' , ||u|| L 2 (Q/) < C||e(u) 



\LHQ f )> 



(iv) For any u £ H 1 (Q f,M. N ) with both Jg u(x) dx = and x x u(x) dx — 0, one has 
\\u\\ L 2 (Qf) < C||e(u)|| L 2 (Q/) . 

In the last assertion, x x u is the skew-symmetric matrix (xiUj — XjUi)^^, and the condition 
means that u is orthogonal (in L 2 ) to the rigid motions (of the form a + Bx with B skew- 
symmetric) . 

Let us sketch the proof of this proposition. First of all, the proof of point (i) is standard (by 
integration along vertical lines starting from Q') and it is well known that the constant C, in this 
case, only depends on M (/ could then be any l.s.c. function below M). In the same way, the 
proof of (iii) is significantly simpler than the proof of (iv) (note however that it does require that 
/ is Lipschitz and C will depend on both M and L), and we will not discuss it. (See [H] for a 
detailed proof, in dimension two). 

To prove (ii) one first establishes the following inequality: there exists Cq depending only on 
M such that for any / and any u £ H 1 (Qf), one has 

(B.l) f u{xfdx < C ( / u(x) 2 dx+ [ |Vw(x)| 2 ) dx , 

JQ f \Ja JQ f J 

where A denotes the set Q' x (0, 1) (the important fact here being that A is an open set that 
belongs to all the domains Qf, for all admissible /). The proof of ljB.l|l . again, is standard. It 
relies on integration along vertical lines starting from the base Q' and on the obvious fact that for 
any b £ [1, M] and any v £ C^O, 6), 

rv /■! pb 

(B.2) / v{t) 2 dt < 2M v{t) 2 dt + 2M 2 v'(t) 2 dt. 

Jo Jo Jo 

Now, if (ii) is not true, it means that there exists functions /„ and u n with 1 < /„ < M, f n 
L-Lipschitz, u n £ H 1 (Qf n ), Jq^ u n dx — and 

ll^nlli^Q^) > n\\Wu n \\ L 2 {Qf j 

for any n. 

Without loss of generality we may renormalize u n so that J A u 2 l dx — 1 . Then, by IjB.ip we find 

1 ^/~Co 
\\^ u n\\^(Qf n ) < -|l u ™IU 2 (Q/„) + < — — ( 1 + I|Vu„||l2 (Q/ j) . 

If we extend both u n and Vu ra with the value outside of Qf n , this inequality shows that V«„ 
goes to zero strongly in L 2 (Qm',^ n ) (Qm = Q x (0, M)), and, up to a subsequence, that there 
exists u £ L 2 (Qm) such that u n — u weakly in L 2 (Qm)- 
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On the other hand, by Rellich's theorem, u n — > u strongly in L 2 (A) and since V« = in A and 
J A u 2 — 1, u is the constant ±1/ yplj". 

Now, since the functions /„ are uniformly equibounded and equicontinuous, up to a further 
subsequence, we may assume also that f n converges to some / uniformly. It is now easy to check 
that u 6 H 1 (Q f), u = outside of Qf, and Vu = (the limit of Vw n ) in Qf so that u is a constant 
in Qf. We deduce that u = (l/^\A\)xQ f ■ Now, for each n, one had u n dx = Jq m u„ dx = 0, 

hence in the limit Jg^ udx = 7^ a contradiction. Hence (ii) must be true. 

We observe here that (ii) holds in fact as long as / belongs to a fixed set of functions which 
is compact in C°(Q', [1, +00)) (the constant C depending only on this compact set). The case of 
L-Lipschitz functions uniformly bounded by the constant M is a particular case. On the other 
hand, for the Poincare-Korn inequalities (iii) and (iv), the fact that the functions / are uniformly 
Lipschitz seems to be essential, as we now show. 

Let us now prove (iv). It is enough to show that the vectorial version of IjB.ll) holds, that is, 



(B.3) 




This will be shown, again, by integration along lines and using (|B.2I) . however, this time, it is not 
sufficient to consider only vertical lines starting from Q' . Indeed, one has for any smooth vectorial 
field u € C-iQf) that 

d(u(x + sO-0/ds = (e(u)(x + sQQ-Z 

for any x £ Qt, £ S S^ -1 , and sd such that x + s£ € Qf. Hence, integration along vertical 
lines will control the component un of u = (ui, . . . , un)- To control the other components, one 
needs to integrate along lines in at least N — 1 other independent directions (as is done in [oj). 
Given p > max{L, 2M } let us consider, for any i = 1, . . . , N — 1, the vectors 

g= J— 7 (0,...,0,±l,0,...,p) 
V 1 +P 

where ±1 appears at the zth position. Given i < N — 1, if we considers the lines starting from 
Q' in the direction we see that they "see" all points x € Qf with Xi > 1/2. On the other 
hand, the lines starting from Q' in the direction £~ "see" all the points with Xi < 1/2. Integrating 
along these lines and using i|B.2() . one controls the L 2 -norms on one half of the domain Qf of 
(ui + piijy) /V 1 + P 2 and, on the other half, of {—Ui + pu^) /y/l + p 2 - Together with the control 
of Jq u 2 n dx obtained previously, this shows that one can control Jq u 2 dx with the right-hand 
side of (|B . 3|) . Repeating this argument for all i, we find that l|B.3|) holds, with now a constant 
that depends on M and L, through p. 

We deduce, exactly as before, that (iv) holds. □ 
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